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Abstract 

Active motion is a concept in complex systems theory and was success- 
fully applied to various problems in nonlinear dynamics. Explicit studies 
for gravitational potentials were missing so far. We interpret the Fried- 
mann equations with cosmological constant as a dynamical system, which 
can be made 'active' in a straightforward way. These active Friedmann 
equations lead to a cyclic universe, which is shown numerically. 

1 Active Motion 

To account for self-driven motion as observed in biological systems, an additional 
degree of freedom, called the "internal energy" e, was introduced in the context 
of complex systems theory ([T], [5]). First interpretations of this formalism 
were given with respect to animal movement and complex motion of particles 
in general. It describes particles who can convert their internal energy e into 
mechanical energy and thus exhibit complex movement. Various applications 
of this model were studied extensively, e.g. in swarm theory [3]. Dynamics of a 
particle with position q and momentum p undergoing active motion as discussed 
in ^ are 

Qi =Pi, Pi = (1 - e) , (1) 

where the evolution of the internal energy e is given by the equation 

e = ci - c 2 e - c 3 eU(q). (2) 

The Cj remain constants to be fixed. The first integral of the active dynamical 
system reads 

f + Ufa) - f dt' e(t')U( qi (t')) = C= const, (3) 
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which can be viewed as a generalization of the classical energy equation. In the 
context of active motion the particle is not simply driven to the minimum of 
the potential U(q), its internal energy allows it to exhibit self-driven movement. 
Depending on the structure of the potential, the systems shows various complex 
motion patterns, which can be studied analytically by bifurcation theory and 
numerically by just simulating the evolution of q, p and e in time. So far, explicit 
investigations were done only for harmonic potentials (see [I] for a detailed study 
of the nonlinear dynamics) . One can naturally ask for active motion of particles 
driven by gravitational interactions. We will now give a formulation of this 
problem within the framework of Friedmann equations. 



2 Active Friedmann Equations 



The Friedmann equations for the dimensionless normalized cosmological scale 
parameter a(t), describing the spatial evolution of a flat, homogeneous and 
isotropic universe with cosmological constant A read 
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where fl m ,o and Q r ,o ar e the density parameters of matter and radiation at 
present time t = and Cl\ = ^rri is the density parameter of the cosmological 

constant. <z(0) = 1 per definition and Ho — d(0) denotes the current Hubble 
parameter. If we introduce the following potential 

U{a) = — — \—t+ +n A a , (6 
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then the analogy with a classical dynamical system of a single particle with 
coordinate a(t) and energy C = is evident. According to general relativity the 
constant C is related to the spatial curvature of the universe which, however, is 
taken to be zero due to recent observations. The density parameter of A is fixed 
by the Friedmann equation ^ which at present time t = gives the relation 

f2 r ,o + ^m.o + Ha = 1- (7) 

We now want to propose a new phenomenological model for the development 
of the universe which is based on the active generalization of the conventional 
Friedmann equations. The active Friedmann equations are given by 



h = -Ht ^i + ^-Ova (1-e), (8) 
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C i- C2 e + c 3 e^(^ + ^+fi Aa 2 ), (10) 



where ci, C2, C3 are arbitrary but fixed real constants. For notational convenience 
we have chosen HqC as the corresponding constant for the first order integral 
in ([9]). Differing from the conventional case we allow for a non- vanishing C in 
the active context, which has to be determined yet. We first have to assign spe- 
cific values to the density parameters of matter/radiation and the cosmological 
constant. Notice that we can't set £l\ = 0.7, since this results from equation 
([7]), which looks different in the active formulation, namely 

fl\ = 1 — Q r fi — Q m ,o — C, (11) 
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which is derived from equation by setting t = 0. The constant C can 
be fixed by demanding that the current acceleration a(0) of the universe in the 
active scheme should agree with the acceleration value given by the conventional 
Fricdmann equations (thus accounting for the experimental result a(0) > 0). For 
any fixed eo := e(0) one finds 

C = (ln m . + 2fl rfi - . (12) 

1 - e V 2 / 

Accordingly, 

n A = — *— (l - n m , (i + ?) - nr,o(l + eo)) , (13) 

1 — e V 2 J 

so that J7a is now completely determined by the density parameters of mat- 
ter/radiation and the initial condition for e, which is freely chosen. Depending 
on this initial condition £l\ may also be brought to vanish. 



3 Discussion 



Figure 1 shows a simulation example of the active dynamics of the normalized 
cosmological scale parameter a(t) . We observe oscillatory behavior with signifi- 
cantly smaller amplitudes and shorter periods in the past. The scale parameter 
a(t) is oscillating around a certain equilibrium point, which can be calculated 
analytically. The active Fricdmann equations can be rewritten as follows: 



P 

dU(a) 



(1-e) (14) 



da 

e = c\ — c%e — c^eU (a) . 
Equilibrium points are found by searching values (a,p,e), for which the right 



hand side of ( 14 1 vanishes. Equilibrium points with e — 1 are unstable, the case 



1 leads to the conditions 

p = 0, 



= fl^a 4 — fl r fi, (15) 
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For the specific choice of parameters and initial conditions made for the sim- 
ulation shown in Figure 1, the numerical values of (a,p,e) can be calculated 
directly. The quartic equation for a has only one positive, real-valued solution, 
namely 5 = 1.01, resulting in (a,p,e) — (1.01,0,4.08). Figure 1 shows that this 
value of the scale parameter marks the center point of the oscillations. The 
stability analysis of the above equilibrium point (a,p,e) can be made by lin- 



earizing the system ( 14 1 and calculating the eigenvalues of the Jacobian matrix. 
Our calculation shows the presence of a purely imaginary pair of eigenvalues 
Ai.2 = Hu) and we are near a Hopf bifurcation point (for bifurcation theory, 
see [5]). In our case w = 1.19, so that the period of oscillations is given by 
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Figure 1: Numerical result for the evolution of a(t), with c\ = 5, C2 = 1, 
c 3 = -1, e(0) = 100, H = 1, f2 ro>0 = 0.29 and fi r ,0 = 0.01. (These are rough 
representative values for the density parameters, consult [5 for a recent overwiev 
of cosmological parameters.) This choice of parameters leads to H,\ — 0.15. 

T = = 5.28 in units of the Hubble time H^ 1 , which agrees well for large 
times with our numerical simulation. 



Summarizing, we can say that interpreting the Friedmann equations as a 
nonlinear dynamical system and handling it within the framework of active mo- 
tion results in oscillatory solutions for the scale parameter. Hence we succeeded 
in constructing a model for a cyclic universe giving rise to a sequence of expan- 
sions and contractions without any singularity, yet accounting for the observed 
spatial flatness and the current accelerated expansion. Besides this particular 
example, the active Friedmann equations allow for a lot of other possible sce- 
narios which deserve further analysis. Whether the active universe represents a 
reasonable model for cosmology will depend on finding a physically meaningful 
interpretation of the additional parameter e. Attempts can be made to identify 
its effect with corrections arising e.g. from cosmological theories with extra di- 
mensions, which will be the task of future investigations. 
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